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Low-temperature specific heat of tlie monolayer high- Tc superconductor Bi2Sr2-a:La3;Cu06+« has 
been measured close to the optimal doping point (a; ~ 0.4) in different magnetic fields. The identi- 
fication of both a term in zero field and a ^/H dependence of the specific heat in fields is shown 
to follow the theoretical prediction for d-wave pairing, which enables us to extract the slope of the 
superconducting gap in the vicinity of the nodes {va, which is proportional to the superconducting 
gap Ao at the antinodes according to the standard d^2_y2 gap function). The «a or Ao (~ 12 meV) 
determined from this bulk measurement shows close agreement with that obtained from spectroscopy 
or tunneling measurements, which confirms the simple d-wave form of the superconducting gap. 

PACS numbers: 74.72.Gh, 74.20.Rp, 74.25. Bt 



I. INTRODUCTION 

In the study of high-Tc cuprate superconductors, low- 
temperature specific heat (LTSH) has proven to be a 
helpful tool in identifying the d-wave symmetry of the 
superconducting gapi^ In zero magnetic field, the quasi- 
particle density of states of a d-wave superconductor 
shows a linear energy dependence owing to the pres- 
ence of line nodes, which is expected to give rise to a 
temperature-dependent electronic specific heat. In 
magnetic field H, an energy shift (Doppler shift) to the 
nodal quasiparticle spectrum becomes important, and 
the electronic specific heat of a d-wave superconductor 
is predicted to have a characteristic y/HT dependence."^ 
In LTSH experiments on high- Tc cuprates, although it 
still seems to be somewhat controversial in the identifica- 
tion of the electronic specific heat in zero field (partly 



owing to its small magnitude) ^i^i^ the -s/^T-dependent 
specific heat in fields has been widely recognized in both 
YBasCugOT-a (YBCO) and La2-^Sr^Cu04 (LSCO)^i^ 
hence providing the bulk evidence for d-wave pairing in 
these prototypical high-Tc compounds. 

Recently, more detailed work showed that a key pa- 
rameter of the d-wave superconducting gap, i.e., the gap 
slope in the vicinity of the nodes , also can be quanti- 
tatively determined from LTSHj^^i This makes LTSH an 
effective measure of the superconducting gap maximum 
at the antinodes, Aq, which is proportional to the ua ac- 
cording to the standard (1^2 _y2 gap form A = Ao cos(2(/)), 
as exemplified by experiments on LSCO across a wide 
doping range. ^ The determination of the superconducting 
gap from LTSH has the merits that it is sensitive to the 
gap structure near the nodes by probing the low-energy 



quasiparticle excitations and it refiects the bulk property 
of the sample. A comparison of it with more traditional 
gap measurements such as spectroscopy or tunneling may 
give valuable insights into the superconducting state. 

In light of the above, we present in this paper the LTSH 
experiment on a monolayer high-Tc cuprate supercon- 
ductor Bi2Sr2_2:LaxCu06+5 (La-Bi2201).Up to now, the 
field-dependent LTSH study to identify the d-wave pair- 
ing effect has seemed to be lacking for this material.!^ 
On a single-crystal sample with a doping level close to the 
optimal {x ^ 0.4), we have resolved a T^ electronic spe- 
cific heat in zero field and a \/i?T-dcpcndcnt electronic 
specific heat in fields, both in conformity with the d-wave 
pairing theory. This has then allowed us to quantitatively 
determine the v/^ and a superconducting gap amplitude 
Ao of ^12 meV for the sample. The wa or Aq is found 
to show reasonable agreement with that determined by 
angle-resolved photoemission spectroscopy (ARPES) or 
scanning tunneling microscopy (STM) for La-Bi2201 at 
the optimal doping level4"— 



II. EXPERIMENT 

The La-Bi2201 single crystal with x ~ 0.4 and a mass 
of 10.2 mg was grown by the traveling solvent floating- 
zone method as reported elsewherei>iS The Tc of the sam- 
ple is 28 K, defined by the onset of the magnetization in 
an ac susceptibility measurement The LTSH measure- 
ment was performed using a thermal relaxation method, 
as described in detail previously, It was done in an Ox- 
ford Maglab cryogenic system in the temperature range 
2.5 — 10 K and in magnetic fields up to 12 T [H \\ c 



2 



250 




qI . , . , . , . , , ,_ 

2 4 6 8 10 



T(K) 

FIG. 1. (Color online) Specific heat of the sample in zero 
magnetic field, plotted as C/T vs T (squares). The solid line 
is the fit to C = l3oT^ + rjoT^. The inset shows the specific 
heat below 6.5 K in zero field and in 12 T, plotted as C/T vs 



axis). The heat capacity of the addenda was measured 
separately and was ^30% of that of the sample in the 
measured temperature range. The chip thermometer has 
been calibrated in fields and the heat capacity of the ad- 
denda is mostly due to phonons and hence shows nearly 
no field dependence.— 



III. DATA AND ANALYSIS 

Figure [T] shows the specific heat of the sample in zero 
field. We found it could be well fitted by the formula 
C = (3oT^ + r]oT^ with (3o and ijo representing the coef- 
ficient of the phonon term and term respectively. 
The fit is shown by a solid line in Fig. [T] and Table |T] 
lists the parameters /3o and rjQ. This indicates that in 
zero field the specific heat comes mainly from phonon 
excitations. Note that the /3o for the present sample is 
roughly five or eight times larger than that reported for 
optimally doped YBCO or LSCO, respectively ; -^^i-^'' re- 
fiecting a substantially greater phonon specific heat for 
La-Bi2201. From (3o the Debye temperature O^i is calcu- 
lated as 226 K. A residual linear specific heat, 7(0)T, is 
usually observed for YBCO and LSCO in zero field^'^'i 
while for our La-Bi2201 crystal, the above fit suggests 
that such a term should be close to zero. This can be 
further illustrated by the inset of Fig. [T] In plot of C/T 
vs T^, it is seen that the 0-T data show a nearly zero 
intercept on the C/T axis at T = K, indicating a neg- 
ligible 7(0)T term. It is worth noting that a nearly zero 
7(0) has also been reported previously for La-free Bi2201 
and BiaSraCaCuzOg+s (Bi2212) single crystalsJ^ The 
reason why the 7(0)r is shown for YBCO and LSCO 
while not for most Bi-based cuprates at the optimal dop- 
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FIG. 2. (Color online) Global fit (solid lines) to the specific 
heat in fields up to 12 T (symbols), plotted as C/T vs All 
data are fit simultaneously to Eqs. ([1} and ([2]) with /3 and rj 
kept the same for all H . For clarity, the interval between data 
sets at different H has been enlarged and the in-field data and 
fit are shown as C/T = 20y (H) + ^T^ + riT^ . The inset shows 
the difference between the data and the fit, (Cat — C)/C vs 
T. 

ing level is presently not clear. ^ In the inset of Fig.[T] we 
have also shown the specific heat in a field of 12 T. A 
nearly parallel shift to the zero-field data, though small, 
can be readily seen. This gives a finite intercept on the 
C/T axis for 12-T data and hence implies that there is a 
field-induced linear-T specific heat, j{H)T. 

One may note that no aT^ term expected for a d- 
wave superconductor is contained in the above fitting to 
the zero-field data. This is because we found that if we 
allowed it to be in the fit, we would obtain a negative 
value of the coefficient a. In fact, previous studies have 
shown that it is usually difficult to identify the T^ term 
by fitting the zero-field data alone because this term is 
expected to be rather small and could be easily concealed 
by the large phonon specific heat .^fe' — Using the data in 
different fields together to more accurately determine the 
phonon contribution is required to give a reliable evalu- 
ation of the aT^ termjS, Therefore, to correctly identify 
the T^ term in zero field and determine the linear-T spe- 
cific heat in fields, we have performed a global fit to all 
the LTSH data,^^ in which the zero-field data set is fit to 

C(r, 0) = aT^ + /3T^ + riT^ (1) 

and the in-field data sets {^J,QH > 0.5 T) are fit to 

C(T,H) ^j{H)T + PT^ + riT^. (2) 
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TABLE I. Fit parameters in the present LTSH study. /3o 

rio [t]), a, and 'y{H) are in units of mj moP^ K^" with n = 4, 

6, 3, and 2, respectively. 



Fit 


Parameters 


H = 


Po = f .86 ± O.Of , 770 = 1.14 X 10"^ 


Global fit 


/3 = 1.84 ± 0.01, n = 8.93 x 10"'', 




a = 0.10 ± 0.04, -f{H) = (0.90 ± 0.06)/ff 



In this global fit, all data sets are fit simultaneously to 
Eqs. (H)) and ^ with a weighted-least squares criterion 
[minimizing X]v "c"'^' )^]'^ The phonon coefRcients /3 
and rj are kept as the same for all data while no assump- 
tions on the magnitude of a or on the field-dependent 
form of j{H) are made. Figure [2] presents the global-fit 
result, where, for clarity, the interval between data sets 
at different H has been enlarged 20 times according to 
the fit. One can see the data are reproduced well by solid 
fit lines, showing the reasonable quality of the fit. The 
inset of Fig. [2] shows the difference between the data and 
the fit, from which the total rms deviation of the data 
from the fit was calculated to be 1.1%. 

The parameters yielded from this fit are listed in Ta- 
ble m with 90% statistical confidence levels. Unlike the 
single fit to zero-field data, the a determined from the 
global fit is positive, supporting the existence of a 
term in zero field. It is noted that the size of the a 
is comparable to that determined for optimally doped 
YBCO or LSCO J^iii Figure [3] plots the -f{H) at each 
field. Plotted together is a solid line showing the fit to 
"f{H) = A\fH based on the theoretical prediction for 
a d-wave superconductor.'^ It is seen that the solid line 
gives a reasonable description to the "f{H), suggesting 
that the field-induced specific heat essentially follows the 
predicted \/HT dependence, in line with the observation 
in YBCO and LSCO.— i^i^ This, alongside the presence of 
an aT^ term in zero field, provides evidence for the dom- 
inance of d-wave symmetry in the superconducting gap 
in optimally doped La-Bi2201. Note that although the d- 
wave pairing has been probed for La-Bi2201 by other ex- 
perimental techniques such as ARPES,— present LTSH 
findings confirm that it is a bulk property of the sample. 

A more quantitative analysis can be made. Within a d- 
wave framework, the magnitude of the aT^ term depends 
on the density of states near the gap nodes, and hence 
correlates with the nodal gap slope, wa, which character- 
izes how steeply the gap opens around the nodes. Specif- 
ically, the coefficient a can be expressed as 

18C(3) fc| nKnoi 1 
a = TT—i ' (3) 

■K h d VFVA 

where C(3) — 1-2, fcs is Boltzmann's constant, h is 
Planck's constant, n is the number of Cu02 planes per 
unit cell, d is the unit-cell size along the c axis, Vmoi is 
the mole volume of the unit cell, and vp is the Fermi ve- 
locity at the nodesii'^'2^ Similarly, as the field-induced 
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FIG. 3. (Color online) Coefficient of the field-induced linear- 
T specific heat, 7(-H"), obtained from the global fit (circles). 
Error bars are at a statistical 90% confidence level. The fine 
is the fit to 7(_H") = {A = 0.90 mJ mol"! K'^ T'^-^) 

for d-wave pairing. 

specific heat A-JTlT of a d-wave superconductor comes 
from the Doppler shift of the quasiparticle states (owing 
to the superfiuid flow around the vortices) in the vicinity 
of the nodeSf^ its magnitude is also sensitive to the da. 
The relation between the prefactor A and the f a is given 

by 

3h \ <^o d va' ^ ' 

where $o is magnetic flux quantum and a — 0.465 for a 
triangular vortex lattice»ii2iiS^ Equations ([3]) and (HI) ex- 
plicitly show that with the a or A determined in the 
experiment, one can in principle extract the va from 
LTSH with no adjustable parameters, which is quite 
helpful because, according to the standard d-wave form 
A = Aocos(2(/)), one can further obtain a more familiar 
quantity, i.e., the superconducting gap maximum Aq at 
the antinodes, via 2Ao — hkpVA, where kp is the Fermi 
wave vector along the nodal direction (0, 0) — (tt, tt). 

For La-Bi2201, n = 2, d = 24.4 A, and Kioi = 
106.8 cm^. A recent ARPES experiment on optimally 
doped La-Bi2201 showed that along the nodal direction 
VP ~ 2.7 X 10^ m s-\ and kp ~ 0.74 A^^.^s ygj^g 
these parameters and the a or ^ as shown in Table [H 
we obtain Aq ~ 13 ± 5 nieV or ~ 10.4 ± 0.7 meV ac- 
cording to Eq. ([3]) or ([1]), respectively. It is seen that 
both determinations of the Aq, which are largely inde- 
pendent, show a fairly good consistency, although the 
former has a larger uncertainty owing to the larger un- 
certainty in determining the a. This further validates 
the analysis of the LTSH data within the d-wave pairing 
theory. From Aq ~ 10.4 ± 0.7 meV, which is determined 
from the field dependence of the specific heat and hence 
has a higher degree of accuracy, one can have a gap-to- 
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Tc ratio Aq /ksTc ~ 4.3 ± 0.3 for our sample, which is 
larger than the weak-coupling d-wave BCS prediction^! 
Ao/fcsTc — 2.14 and comparable to the values found in 
other optimally doped high-Tc cuprates such as YBCO 
(- 4.3 in Ref.^), LSCO (~ 4.9 in Ref. fl), and Bi2212 
(- 4.5 in Ref. 123. 



IV. DISCUSSION 

Focusing on La-Bi2201 at the optimal doping level 
(x ~ 0.4), it is instructive to compare the Aq determined 
here to what has been measured by other experimen- 
tal techniques such as ARPES and STM. In an early 
ARPES experiment^ an antinodal superconducting gap 
of 10±2 meV was observed, which shows good agreement 
with the LTSH result. This agreement confirms that the 
superconducting gap essentially follows the simple d-wave 
form because in ARPES the Aq is measured directly at 
the antinodes while in LTSH it is derived from the nodal 
gap slope VA with the assumption A = Aocos(20). In 
several recent ARPES studies, it was found that from 
nodes to antinodes the superconducting gap can indeed 
be well fitted b y t he simple d-wave form with Aq of 
13.5 meV (Refs."|l3 and[lll) or 15.5 meVJ^ In a recent 
STM experiment an averaged superconducting gap of 
11.4 ± 4 meV was determined from spatial conductance 
maps. We can see these reports on Aq agree quantita- 
tively with the LTSH result as well. In Fig. |4] we have 
plotted the Aq for near optimal-doped La-Bi2201 inferred 
from LTSH [according to Eq. (j?])] and aforementioned 
ARPES and STM studies, S.ii2,"i^ where the doping level 
of the sample (p ~ 0.16) is determined by the empirical 
relation^liS^ between p and the reported La content. 

Unlike the above-compared ARPES experi- 
ments,— i^^rJ^ Kondo et al. reported an antinodal 
gap as large as ^40 meV below Tc in optimally doped 
(Pb,La)-Bi2201 and ascribed it to the presence of a pseu- 
dogap with a nonsuperconducting origin.^ To compare 
LTSH with such an ARPES experiment where a different 
energy scale other than the superconducting gap seemed 
present at the antinodesj^ one needs to focus directly 
on the VA near the nodes. In fact, in the work of Kondo 
et al, it was shown that an extrapolation of the gap 
measured around the nodes under the standard d-wave 
form would yield a gap amplitude of ~15 meV at the 
antinodeSf? suggesting that the ua they probed is also 
comparable to that obtained in LTSH. This further helps 
reveal that, whether a large pseudogap was detected 
or not at the antinodes below Tc, in recent ARPES 
measurements^ri^ the superconducting gap determined 
near to the nodes exhibited consistent behavior, i.e., a 
quite similar nodal gap slope. 

The above comparison shows that LTSH provides a 
reliable determination of the superconducting gap in La- 
Bi2201 at the optimal doping level, mimicking our pre- 
vious observation in LSCO at various doping levels. As 
pointed out earlier, in recent work we have also deter- 
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FIG. 4. (Color online) Comparison of the superconducting 
gap Ao in near optimally doped La-Bi2201 obtained from the 
present LTSH study [according to Eq. Q, soli d sq uare] and 
from some ARPES (solid circles) (RefsHl, [MSI) or STM 
(solid triangle) (Ref. [ll|) studies. The Aq from LTSH is in- 
ferred from the nodal gap slope va under the standard d-wave 
form A = Ao cos(2(/)). The Ao from ARPES was measured at 
the antinodes in Ref. ^ or determined with the observation 
of a standard d-wave superconducting gap in Refs. 40- 
To gain a better perspective, the comparison of the Ao in 
LSCO obtained from LTSH (crossed squares) (Refs.@ andQ), 
ARPES (crossed circles) (Ref. [H), and STS (crossed trian- 
gles) (Ref. I30I ) is also plotted in the doping range from slight 
overdoping to underdoping. The dashed line is a guide to the 
eye. See the text for details. 



mined the superconducting gap in LSCO across a wide 
doping range by probing the field dependence of the 
LTSH.^'^ In the overdoped regime, the obtained Aq was 
found to show a quantitative agreement with that derived 
from other experimental techniques such as scanning tun- 
neling spectroscopy (STS)j^ In the doping range from 
slight overdoping to underdoping. Fig. 3] shows that in 
LSCO the Aq obtained from LTSH is also well confirmed 
by STS or recent ARPES measurements , ^^i Hence this 
illustrates the applicability of extracting the supercon- 
ducting gap from LTSH across nearly the whole super- 
conducting phase diagram of LSCO. Note that it also in- 
dicates a growing of the va below optimal doping toward 
the underdoped regimej^*- and that the superconducting 
gap basically retains the simple d-wave form with dimin- 
ishing p (at least down to p ^ 0.11), which is of impor- 
tance particularly in the context of a recent debate on 
the gap property of the underdoped high- Tc cuprates 4^ 
In view of these findings in LSCO, we believe it would be 
appealing in the future to explore the LTSH of La-Bi2201 
at more doping levels, with the expectation of acquiring 
the doping dependence of the superconducting gap in this 
archetypal cuprate from a bulk, thermodynamic method. 
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V. CONCLUSION 



In summary, we have measured the LTSH in different 
magnetic fields in near optimally doped La-Bi2201. By 
performing a global analysis of the data, we have identi- 
fied a electronic specific heat for H = and a ^/HT 
electronic specific heat for iJ 7^ 0, as suggested for a d- 
wave superconductor. This provides the bulk evidence 
for d-wave pairing in optimally doped La-Bi2201, joining 
previous findings in YBCO and LSCOiii^ii Moreover, we 
have quantitatively determined the slope of the super- 
conducting gap Vi^ near the nodes and its amplitude Aq 



(~ 12 meV) at the antinodes according to the standard 
dx2_y2 form, which show close agreement with the re- 
sults of ARPES and STM.^-i^ As shown in LSCO,-^-! the 
present study further demonstrates the virtue of LTSH 
as a bulk method to probe the superconducting gap near 
the nodes in high- Tc cuprates. 

ACKNOWLEDGMENTS 

Technical assistance from Jing Wang is gratefully ac- 
knowledged. This work is supported by NSFC, CAS 
(Project ITSNEM), and MOST of China (973 Program, 
Nos. 2006CB601000 and 2011CBA00100). 



* |yue . wangQp ku . edu . cn| 
t hhwen@nju.edu.cn 

1 N. E. Hussey, Adv. Phys. 51, 1685 (2002). 

^ R. A. Fisher, J. E. Gordon, and N. E. Phillips, in Hand- 
book of High- Temperature Superconductivity: Theory and 
Experiment, edited by J. R. Schrieffer and J. S. Brooks 
(Springer, New York, 2007) Chap. 9. 

^ G. E. Volovik, JETP Lett. 58, 469 (1993). 

* J. W. Loram, J. Luo, J. R. Cooper, W. Y. Liang, and J. L. 
Tallon, J. Phys. Chem. Solids 62, 59 (2001). 

^ Y. Wang, B. Revaz, A. Erb, and A. Junod, Phys. Rev. B 

63, 094508 (2001). 
^ H.-H. Wen, L. Shan, X. G. Wen, Y. Wang, H. Gao, Z.-Y. 

Liu, F. Zhou, J. W. Xiong, and W. X. Ti, Phys. Rev. B 

72, 134507 (2005). 
^ Y. Wang, J. Yan, L. Shan, H.-H. Wen, Y. Tanabe, 

T. Adachi, and Y. Koike, Phys. Rev. B 76, 064512 (2007). 
^ J. M. Harris, P. J. White, Z.-X. Shen, H. Ikeda, 

R. Yoshizaki, H. Eisaki, S. Uchida, W. D. Si, J. W. Xiong, 

Z.-X. Zhao, and D. S. Dessau, Phys. Rev. Lett. 79, 143 

(1997). 

^ T. Kondo, T. Takeuchi, A. Kaminski, S. Tsuda, and 
S. Shin, Phys. Rev. Lett. 98, 267004 (2007). 

'° J. Wei, Y. Zhang, H. W. Ou, B. P. Xie, D. W. Shen, J. F. 
Zhao, L. X. Yang, M. Arita, K. Shimada, H. Namatame, 
M. Taniguchi, Y. Yoshida, H. Eisaki, and D. L. Feng, Phys. 
Rev. Lett. 101, 097005 (2008). 

" J.-H. Ma, Z.-H. Pan, F. C. Niestemski, M. Neupane, Y.- 
M. Xu, P. Richard, K. Nakayama, T. Sato, T. Takahashi, 
H.-Q. Luo, L. Fang, H.-H. Wen, Z. Wang, H. Ding, and 
V. Madhavan, Phys. Rev. Lett. 101, 207002 (2008). 

^2 J. Q. Meng, W. T. Zhang, G. D. Liu, L. Zhao, H. Y. Liu, 
X. W. Jia, W. Lu, X. L. Dong, G. L. Wang, H. B. Zhang, 
Y. Zhou, Y. Zhu, X. Y. Wang, Z.-X. Zhao, Z. Y. Xu, C. T. 
Chen, and X. J. Zhou, Phys. Rev. B 79, 024514 (2009). 
B. Liang and C. T. Lin, J. Cryst. Growth 267, 510 (2004). 
The actual La content of the sample is inferred to be 0.36, 
based on the wavelength dispersive spectroscopy (WDS) 
measurements on crystals grown under the same condition 
(Ref. IT^ . With this, the sample's hole doping level p is 
estimated to be ~0.165 according to the empirical relation 
between x and p proposed in the literature (Refs. [IJ and 
[2^ ). This slight deviation of the p from the optimal doping 
point (0.16) may explain the Tc of the sample lower than 



T^'"^'' ~ 32 K. 

^5 H.-H. Wen, Z.-Y. Liu, F. Zhou, J. W. Xiong, W. X. Ti, 
T. Xiang, S. Komiya, X. F. Sun, and Y. Ando, Phys. Rev. 
B 70, 214505 (2004). 

K. A. Moler, D. L. Sisson, J. S. Urbach, M. R. Beasley, 
A. Kapitulnik, D. J. Baar, R. Liang, and W. N. Hardy, 
Phys. Rev. B 55, 3954 (1997). 
" R. A. Fisher, N. E. Phillips, A. Schilling, B. Bufleteau, 
R. Calemczuk, T. E. Hargreaves, C. Marcenat, K. W. Den- 
nis, R. W. McCallum, and A. S. O'Connor, Phys. Rev. B 
61, 1473 (2000). 

A. Chakraborty, A. J. Epstein, D. L. Cox, E. M. McCarron, 

and W. E. Farneth, Phys. Rev. B 39, 12267 (1989). 

J. S. Urbach, D. B. Mitzi, A. Kapitulnik, J. Y. T. Wei, and 

D. E. Morris, Phys. Rev. B 39, 12391 (1989). 

As shown in Fig. [2] there is no discernible Schottky-like 

anomaly in all LTSH data sets at different H , which has 

helped simplify the data analysis. 

M. Chiao, R. W. Hill, C. Lupien, L. Taillefer, P. Lambert, 
R. Gagnon, and P. Fournier, Phys. Rev. B 62, 3554 (2000). 
I. Vekhter, P. J. Hirschfeld, and E. J. Nicol, Phys. Rev. B 
64, 064513 (2001). 

M. Hashimoto, T. Yoshida, H. Yagi, M. Takizawa, A. Fuji- 
mori, M. Kubota, K. Ono, K. Tanaka, D. H. Lu, Z.-X. 
Shen, S. Ono, and Y. Ando, Phys. Rev. B 77, 094516 
(2008). 

^'^ H. Won and K. Maki, Phys. Rev. B 49, 1397 (1994). 
K. Nakayama, T. Sato, K. Terashima, H. Matsui, T. Taka- 
hashi, M. Kubota, K. Ono, T. Nishizaki, Y. Takahashi, 
and N. Kobayashi, Phys. Rev. B 75, 014513 (2007). 

^® H. Ding, M. R. Norman, J. C. Campuzano, M. Randeria, 
A. F. Bellman, T. Yokoya, T. Takahashi, T. Mochiku, and 
K. Kadowaki, Phys. Rev. B 54, R9678 (1996). 
Y. Ando, Y. Hanaki, S. Ono, T. Murayama, K. Segawa, 
N. Miyamoto, and S. Komiya, Phys. Rev. B 61, R14956 
(2000). 

H.-Q. Luo, P. Cheng, L. Fang, and H.-H. Wen, Supercond. 
Sci. Technol. 21, 125024 (2008). 

M. Shi, J. Chang, S. Pailhes, M. R. Norman, J. C. 
Campuzano, M. Mansson, T. Claesson, O. Tjernberg, 
A. Bendounan, L. Patthey, N. Momono, M. Oda, M. Ido, 
C. Mudry and J. Mesot, Phys. Rev. Lett. 101, 047002 
(2008). 



6 



T. Nakano, N. Momono, M. Oda, and M. Ido, J. Phys. Soc. 

Jpn. 67, 2622 (1998). 

A. J. Millis, Science 314, 1888 (2006). 



